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ABSTRACT 

For a particular case of three-body scattering in 2 dimensions, we demonstrate an- 
alytically that the behaviour of the adiabatic potential is different from that of the 
hyperspherical coupling matrix elements, thereby leading to a phase shift that tends 
to zero, as the energy goes to zero, instead of to a constant. We consider two parti- 
cles interacting with binary repulsive step potentials, one acting as a spectator, and 
solve analytically for the adiabatic eigenvalues and eigenvectors, for all values of p. We 
are thereby able to obtain the leading terms in the long range behaviour of the effec- 
tive potentials, and confirm its inverse logarithmic nature, suggested from a previous 
numerical study. 



Introduction 



In a previous paper the authors show how, starting from hyperspherical harmonic 
expansions, they obtained adiabatic potentials, suitable for the calculations of three- 
body phase shifts at low energies. The calculations were for 3 particles in a plane, 
subject to finite repulsive core interactions. 

The calculations were meant to establish a method which would lead to the eval- 
uation, at low temperature, of a third fugacity coefficient in Statistical Mechanics. 
The latter task was subsequently carried out by Jei Zhen and one of the authorsp|. 
In both investigations, it was important to consider different cases, corresponding to 
the distinct representations of the permutation group and different physical situations, 
with either the 3 particles interacting or simply two of them interacting, with the third 
acting as a spectator. 

Absolutely crucial, in these investigations, is the large-p behaviour of the adiabatic 
potentials. The nature of the long "tail" of the adiabatic potential determines how 
the correspondent eigenphase shift behaves, as the energy tends to zero. Thus, our 
most significant result was that for the 3 most important types of the phase shifts, 
associated with the cases of ^Fi^, °r2g and 5, the adiabatic potentials (the adiabatic 
eigenvalue minus a centrifugal term) behave as l/(p^/np), for large values of p, instead 
of the of the hyperspherical potential matrix elements. This then implies that 
the phase shifts, instead of tending to constants as the energy goes to zero, behave as 
l/(ln g)^, and therefore go to zero! (The variables p and q are, respectively, the hyper 
radius and the reduced wave number.) 

Though, in our old paper, our basic material was numerical, we were able never- 
theless to propose "heuristic" formulae, to characterize the asymptotic behaviour of 
the 3 types of eigenpotentials, of the remodeling that takes place to yield a different 
scattering from the one expected from the solution of a finite number of hyperspherical 
equations. 

In this paper, we present analytical results, where we show that in one of the 
three cases mentioned above, the case 6, we succeeded in calculating analytically the 
adiabatic eigenvectors and eigenvalues for all the values of p, and therefore also in the 
asymptotic region. 

While this calculation involves where only two of the particles interact, 

while the 3rd particle acts as a spectator, it is well to note that in the hyperspherical 
coordinate system the two-body interaction is long ranged (in p) and also that in 
the full hyperspherical calculations of the other cases, we only need, using symmetry 
and enforcing a restriction on the quantum numbers, to take into account the matrix 
element of one of the pair potentials. 

Here, then, our calculations allow us to re-examine our previous results, and confirm 
and extend the asymptotic form (and coefficients) that can be used to characterize the 
long range behaviour of the adiabatic potentials. 
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The KL Hyperspherical Coordinate System 
The Harmonic Basis 



For a system of three equal mass particles in two dimensions, we define the Jacobi 
coordinates 



r/= (fi - f2)/V2 and ^=y'2/3( — fg) , 

which allows us to separate, in the Hamiltonian, the center of mass coordinates from 
those associated with the internal motion. 

Kilpatrick and Larsen0 then introduce hyperspherical coordinates, associated with 
the moment of inertia ellipsoid, of the 3 particles, which allows them to disentangle 
permutations from rotations and obtain harmonics which are pure representations of 
both the permutation and the rotation group. Taking the z axis normal to the plane 
of the masses, we write for the cartesian components of the Jacobi coordinates 

Vx = p(cos^9 cos COS'?/' + sin'i? sin sin-^), 

rjy = p(cos "i? cosy) sin — sin ^9 sin (y9 cos '?/'), 

= p(cos 'd sin ip cos ip — sin-d cos (psinip), 

iy = p(cos'!? sin sin?/; + sin-i? cos cos?/'), (1) 

in terms of the hyper radius p and of the three angles (p and if). 
The harmonics, in their unsymmetrized form, are then 

r;;^^]) = CfQ'^f^{x)e''"<'e'^^ (2) 

where x = sin 2{} and 

ef(x) = (l-xr/2(l + xf/^P„"'^(x) (3) 
P"'^(x) is a Jacobi polynomial, and the normalization constant is 

r^afs _J^^ + l^( n + a + P \ I n + a 



20+/3+1 ) \ a ) \ a j 

The hyper radius p satisfies p^ = rf + ^"^ and the angular components have the ranges 

-l<x<l, -7r/2 < < 7r/2, < ip < 27r . 
Finally we have for the indices the relations 

^ = - i^ax{|z/|, |A|}], a = l|z/ + A|, /3 = l|z/-A|, 

where N is the degree of the harmonic, and A is the inplane angular momentum quan- 
tum number. The indices u and A take on the values — to in steps of 2; all three 
have the same parity and = 0, 1, 2, . . . . 

Linear combinations of the basic harmonics can then be formed]^ to obtain irre- 
ducible bases, adapted to the symmetries of the physical problemsjl], 0|. 
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The Adiabatic Basis 



For our model, the particles interact via a binary step potential 

where the height Vq, and the range, a, are both finite. 

The adiabatic eigenfunctions Bi are then defined as satisfying 

1 9m 1 

+ -^V{p,n)^ = Mp)Bi , (5) 

where V{p, Q) is either the sum of the binary potentials or, simply one of the binary 
potentials, say V^ru), expressed as a function of p and the angles. The index / stands 
for the set of quantum numbers which characterize and index the particular class of 
solutions. A;(p) is the eigenvalue, which upon subtraction of a "centrifugal" type term 
yields the effective potential, of concern to us later on. 

The eigenfunctions may now be used to expand the wavefunctions of the physical 
systems: 

<ii = Y.B,ip,n)<pi,{p), (6) 
I' 

where the amphtudes (pi{p) are the solutions of the coupled equations: 

-Y: ldnB!{n,p)—{B,{n,p)<Pi,{p)) + \i{p)<p,{p) = {2mE/h')<Pi{p). (7) 

The adiabatic eigenfunctions can themselves be expanded in hyperspherical har- 
monics and this is how a large set of them were calculated in the papers quoted earlier. 
The symmetries of the hyperspherical harmonic basis are, of course, reflected in the 
solutions of the adiabatic eigenvectors. For the fully symmetric Hamiltonian, the set 
of solutions divides into nine separate subsets^, each requiring calculations involving 
combinations of matrix elements of only one of the binary potentials, but with restric- 
tions on the quantum numbers of the unsymmetrized harmonics involved. In the case 
of two interacting particles, with a third as a spectator, we flnd an additional four 
subsets. 

The numerical approach was then, for each p, to evaluate a large potential ma- 
trix, with the appropriate harmonic basis, add to this the (diagonal) "centrifugal" 
contribution arising from the angular part of the kinetic energy (the angular part of 
the Laplacian in the Hamiltonian) and diagonalize to obtain the required adiabatic 
eigenvalues. The number of harmonics, needed for numerical convergence, increases 
as a function of p, but it was our fortunate experience to flnd that it was possible 
to evaluate correctly the eigenvalues, that we sought, for values of p large enough 
that the behaviour of A/(p) could be described by asymptotic forms. We were able to 
characterize them, and this gave us the values of Xi{p) for all the larger values of p. 
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Dual Polar Set of Coordinates 



The Harmonic Basis 

In this part of the paper we wish, exclusively, to consider the case of two particles 
interacting together, the third acting as a spectator. As we shall show, we are then 
able to obtain exact adiabatic solutions. 

Our reasoning is as follows. When the third particle does not interact with the 
other two, this must imply that the motion of the pair (1,2), and therefore its angular 
momentum, is unaffected by the motion of the third particle. In a parallel fashion, the 
motion of the third particle, and its angular momentum about the center of mass of 
the particles (1,2), must be a constant as well. If we choose our coordinates carefully, 
the angular behaviour of two of the angles should "factor" out and, for a given p, only 
one variable should be involved in a key differential equation. 

We note that in the KL coordinates, the distances between particles involve two of 
the angles, for example equals p^(l + cos 2i? cos 293) . To get around this, we choose 
an angle to give us the ratio of the length of the 2 Jacobi vectors, and then polar 
coordinates for each of them. Thus, we represent Q by (^i,^2,0), where 77 = pcos0, 
^ = psin0 and rj^ — 77 cos ^1, rjy — 77 sin ^1, — ^^08 02, — Csin^2- The ranges of 
these angles are 

< < 7r/2, < ^1 < 27r, < ^2 < 27r . 

To obtain the harmonics, in a manner which is suitable to also demonstrate the link 
with the KL harmonics, we introduce complex combinations of the Jacobi coordinates, 
i.e. the monomials 



y) 



{Vx - IVy) - Kix - liy) 
Z2 = (77x -'niy)+ - <y) 

Z*2 = {Vx + «%) - «(Cx + <2/) (8) 



It then follows that 



P = l^\ZlZi + Z2Z2) 

and, clearly, zi, zl, Z2 and Z2 each satisfies Laplace's equation, as do the combinations 
ZiZ2-i Z1Z2, zlz2, z^zl and these combinations raised to integer powers. 

Writing p\ — ziz^ and — Z2Z2, we can write as the most general solution arising 
from the monomials zi and Z2'- 

P2 "T Pi 

where Zi, I2 and I are positive integers or zero, and P^'^''^ is a Jacobi polynomial. 
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In terms of the angles, our expression becomes proportional to: 

p'i+'2+2'(cos2 0)'i/2(sinV)'''''P/'''Ucos 2(/))e*^i'ie^^^'^ 

and, finally, in terms of z equal to cos 20, we define our unnormalized harmonic: 

Yl''^\K K = (1 + zf'\l''{l - p/^l'l^il(^) e^«l'le*^2'^ (10) 

where now l\ and I2 can be positive, negative, integers - or zero. (This takes into 
account the other combinations Z\z\^ etc.) The order of the harmonic is N equal to 
+ 1^21 + 2/. 



The Adiabatic Differential Equation 

Writing 

inserting our polar coordinates into the left hand side and changing to our variable 2;, 
we find: 

- 4(1 - .2)_ _ 8.- + + . (12) 

If we now write our adiabatic eigenfunctions as 

Sj^'^(p, Q) = e*^ie*^^ (1 + ^)l'il/2(l - z)l'^l/2 fI'^I'I'^I(p, z) , (13) 
then the functions F will satisfy the equation: 



+ 



{h + l2){h + l2 + 2)+p'V{p,z)]Fl^^'\p,z)=p'Xip)Fl^''\p,z) (14) 



where V(p, z) equals 2m times the potential and in our notation we have dropped 
the absolute value indications. 

When V{p, z) = 0, we can obtain a solution which is analytic between —l<z< +1. 
For A equal to (/i + ^2 + 2/)(/i + /2 + 2/ + 2)/p2 and / a non-negative integer, we find 
that our F is simply pI'^'''^{z), the Jacobi polynomial which appears in our Eq. (10). 
The N that appears in the B of Eq. (13) is the order of the corresponding harmonic. 

For our potential 

V(o z) - I (2W^')V^o -1 < ^ < -1 + l/p' . . 

^^P^^)-\0 -l + l/p2<^<l, ^^^) 

the solutions of this equation which behave reasonably at z equal to —1 and +1 will be 
seen to be proportional to extensions of the Jacobi polynomials to functions with non- 
integer indices, in a relationship similar to that of Legendre polynomials and Legendre 
functions. 
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To motivate and clarify our procedure we first consider the case of li = I2 = 0, with 
and without potential 

When the potential is put to zero and we factor a 4 as well as change the sign, the 
differential equation reads 



d 



d 



l-z')—~2z— + l{l + l) 



dz 



dz 



0. 



(16) 



This is, of course, the Legendre differential equation and, with / a positive or zero 
integer, the well behaved solutions are the Legendre polynomials. 

In the case of our potential, which is zero or a constant (only a function of p) in 
the different ranges of z, we can write our differential equation in a very similar form, 
i.e. as 

{l-z')-"^ ^ ^ 



dz^ '^d-z^"^"^'^ 



0, 



(17) 



where for —l + l/p'^<z<l 



v{v+l)=p'\{p)/A (18) 

and for — 1 < z < — 1 + l/p2 

z/(z/ + l)=p2[A(p)-Fo]/4. (19) 

Denoting the respective values of v as vi and z/2, the corresponding solutions are 
the Legendre function P^^ {z) and the combination 

Py^i-z) = cos(7rz/2) Pi/al-^) " (2A) sin(7rz/2) Quzi^) , 
of the first and second Legendre functions. 

The point is as follows. Whereas Pui{z) is well behaved at z equal to 1, and is 
suitable as a solution for its range in z from — l + l/p^ to 1, both the Pu2{z) and Qu2{z) 
have a logarithmic singularity at z equals —1. The combination that we propose, 
however, is such that the logarithmic terms cancel out and the combination |^ is a well 
behaved solution in the range —1 to —1 + l/p^. 

Expressing these solutions as power series, the first about z = 1, the second about 
z = —1, we obtain 



,Fi(-z/i, z/i + 1; 1; |(1 - z)), for - 1 + l/p^ <z<l 
and 

iFi{-iy2, V2 + 1; 1; \{\ + z)), for - 1 < ;z < -1 + l/p^ 



(20) 



Our overall solutions are then obtained by matching the logarithmic derivative of the 
two solutions (above) at the boundary: at z equal — 1 + l/p^. This then also yields the 
adiabatic eigenvalues. 

It now remains to note that for the cases of l\ and I2 not equal to zero, we can use 
the same procedure. We have, for the two regimes, solutions proportional to 



2Fi(-z/i, i/i + |/i| + I/2I + 1; I/2I + 1; 1(1 - 2)), for - 1 + l/p2 <z<\ 



and 



,Fi(-i/2, V2 + + I/2I + 1; |/i| + 1; 1(1 + z)), for - 1 < z < -1 + l/p2 



(21) 
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For each choice of li and I2 there is an infinite set of values of ui for which the 
logarithmic derivative of the hypergeometric functions can be matched at z equal to 
— 1 + For each such value of i>i, the adiabatic eigenvalue is then given by 

^ (2.. + N + fe| + l)^-l 

When Vo = 0, the adiabatic basis reduces to the hyperspherical harmonic basis 
of Eqn. (10), since the hypergeometric functions reduce to Jacobi polynomials, and 
1/1 = 1/2 = /. So our Sj^''^ is precisely the l^'^'^^l^i; ^2,-2)- 



Comparison of the Adiabatic Eigenvalues 

When the numerical work was done (using the KL basis), lists were made of the appro- 
priate harmonics needed to form the matrices (potential and centrifugal) which, when 
added and diagonalized, yield the adiabatic eigenvalues. We now need to identify these 
eigenvalues and compare them with those obtained by the new method. This is not 
trivial, but an immediate remark can be made. 

First of all, the angular momentum A is a good quantum number, with 

X = h + l2. (23) 

This follows from the fact that li specifies the angular momentum of the 1-2 pair and 
I2 specifies the angular momentum of the third particle relative to the center of mass 
of the first two. Thus their sum defines the total inplane angular momentum. Hence, 
for example, when A = we can have all pairs li and I2 with li = —I2. If h = h = 0, 
this then provides a single eigenvalue for each choice of N = 21, I = 0,1,2, ... . 

Another indicator is wether n is even or odd, which is very significant in the drawing 
up of the lists, associated with the symmetries of the harmonics. Proceeding, then, we 
compare values of the effective potential, defined by 

V{p,N) = X{p)-^-^ i, (24) 

where we subtract from each eigenvalue the value of the centrifugal term that would 
correspond to it, if the binary potential were allowed to go to zero. These have been 
extensively tabulated by Zhenf^. 

Table 1 compares the results in the simplest case, iV = 0, where we demonstrate 
the convergence of the trucated matrix method with the result obtained directly, for a 
sample value of p = 5 and A* = 10. (A* = {li^ /mVoa'^Y/'^) 



N 

^ "max 


\/(5,0) 


110 


0.011754744 


120 


0.011754730 


130 


0.011754670 


140 


0.011754666 


Direct 


0.011754562 



Table 1. Convergence of the matrix method 
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A more extensive set of comparisons is made in Table 2, where selected values of 

the effective potential, obtained from eigenvalues of the truncated matrix, are chosen 
for various values of A^, A and n and compared with the direct results. In all cases, 
except the first, the matrix was truncated at Nmax = 100. 



n 


A 


N 


V{5,N) 


/ 


h 


k 


V(5,Ar) 


E 








0.011754666 








0.011754562 


E 





2 


0.037577818 


1 





0.037577462 








2 


0.000874927 





1 


1 


0.000874911 


E 





4 


0.062609805 


2 





0.062609219 


E 





4 


0.00005971 





2 


2 


0.00005971 





2 


4 


0.00413519 


1 


1 


1 


0.00413512 


E 


1 


1 


0.024168 





1 


0.02416738 


E 


1 


1 


0.00029592 





1 





0.00029591 





1 


3 


0.000024 





2 


1 


0.00002426 





1 


3 


0.00172537 





I 


2 


0.00172529 


E 


1 


3 


0.050462 


1 


1 


0.0504588 


E 


1 


3 


0.00226748 


1 


1 





0.00226737 


E 


2 


2 


0.00000616 





2 





0.00000616 





2 


2 


0.036849 





2 


0.03684737 


E 


2 


4 


0.000088636 


1 


2 





0.000088629 


E 


4 


4 


0.062866 


1 


2 


0.06286247 



Table 2. Some effective potential values in the "delta-bar" class 



Asymptotic Behaviour 

The matching of logarithmic derivatives provides a means of obtaining information 
about the asymptotic behaviour of the eigenvalues, and hence the effective potentials, 
as the hyper-radius, p, gets large. There is however a particular difficulty in finding 
this behaviour. It is that it is not simply a case of looking at the limiting behaviour of 
2^1(0, b; c; e) and 2-^1(0, 6; c; 1 — e) as e — > 0, because the expressions corresponding to 
a and b both depend on p. 

In the simplest case, corresponding to h = I2 = I = 0, we find 



V{p, N) 



Truncated Matrix 



Direct 



4(A^ + l)eo 



1 



(25) 



A + B\np ' 



where 



A 



IoiV2A) , 1 



In 2 



(26) 



2{N + l)V2Ah{V2A) 4(iV+l) 
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and 



1 



2(iV+ 1) 

The Jj's being modified Bessel functions of integer order of the first kind. 
The next simplest case is with only li = 0. Then 



(27) 



A 



1 



(iV + 1) 



2A) 



2V2A/i(V2A) 4 



1 1 1 



p=l 



P 



1 k+m -| 



and 



B 



1 



(2J 



(29) 



2(A^+ 1) ' 

where it should be understood that X)? = 0. 

The case with h ^ introduces considerable complications, particularly to the form 
of A, with higher order Bessel functions occurring, and so is not listed. However the 
expression for B is exactly the same. This confirms the value postulated by Larsen 

It is impressive how well Zhen ^ did working with approximate values for the 
eigenvalues. In her thesis she compares her B^s with the postulated values. If we 
compare her ^'s with the above expressions (see Table 3), we see just how consistent 
her calculations are. 





/ 


k 


A (Zhen) 


A (here) 











2.6064 


2.8293 


2 


1 





0.7581 


0.7764 


4 


2 





0.4146 


0.4159 


1 





1 


1.2381 


1.2897 


3 


1 


1 


0.5493 


0.5511 


5 


2 


1 


0.3356 


0.3327 



Table 3. Comparison of numerical and analytic 
asymptotic leading terms. 

Thus, the conclusions previously obtained by assuming this form of asymptotic be- 
haviour 01 are verified; at least in the "delta-bar" case. 
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Conclusion 



It is now clear that the extensive numerical calculations of Zhen , using the truncated 
matrix approach, provided good estimates of the eigenvalues, the effective potentials, 
and the 2 + 1 phase shifts of the third cluster. The results are consistent for the entire 
range of values of p, taking into consideration the requirement for larger Nmax at larger 
values of p. 

We were also able to demonstrate the all important logarithmic behaviour in the 
asymptotic form of some of the effective potentials. This insures that the corresponding 
phase shifts (dominant at low energies) go to zero, as the wave number goes to zero. 
For the other 2 + 1 phase shifts, characterized by other group classifications of the 
harmonics, we can demonstrate by explicit calculations that both the asymptotic form 
of the effective potentials and the phase shifts go to zero in a stronger manner. 

We would love to obtain similar asymptotic expressions for the effective potentials 
of the fully interacting problem. If we were able to do this, it would simplify enormously 
the cluster calculations, as well as increase its accuracy. 
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